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ABSTRACT 


Rip current spacings are eompared with wave energeties and direetional spreading 
in the Southern Monterey Bay. Southern Monterey Bay affords a unique environment to 
study rip currents owing to their prevalence created by near-normally ineident waves on a 
sandy shoreline. It is hypothesized that rip current spacing is a function of wave 
directional spreading and energy flux, based on the morphodynamic modeling by Reniers 
et al. 2003. A gradient of wave energy flux exists due to headlands and refraction over 
Monterey Canyon. Rip eurrents are shown to occur between cusps in the shoreline, 
allowing cusp spacing to be a surrogate for rip spaeing. Rip current spaeing was inferred 
from beaeh morphology surveys, LIDAR imagery, and Argus cameras, and found to be 
O(150m) at Sand City and O(300m) at Marina, separated by □ 6kni. Measured waves 
during a two month period using wave-rider buoys, show a gradient of aeross-shore 
energy flux between Sand City, Fx = 28000(J/m^), and Marina, Fx = 33000(J/m^) . 
The two sites have the same peak directional spreading of energy value, = 14°, and 

slightly different bulk values for Sand City, cr^„/i.=18°, and Marina, 20°. 

Therefore, the variations in rip current spacing could not be attributed to direetional 
spreading, but appear related to variations in energy flux. 
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I. INTRODUCTION 


A. RIP CURRENTS 

Rip currents are perhaps the most visible feature in the nearshore, in terms of 
scientific interest and public awareness. They are strong, jet-like currents that flow 
seaward across the surf zone in narrow channels, carrying with them sediment up to 
several surf zone widths offshore (Huntley and Short 1992). They are easily 
distinguished by “zones of roiled water and foam,” (Shepard 1936) typically discolored, 
moving offshore and interfering with the incoming waves. Their occurrence is generally 
regulated to intermediate beaches, with reflective and dissipative beaches exhibiting little 
rip activity (Short 1985). Rip currents are an important nearshore process responsible for 
offshore sediment transport, beach erosion and accretion, rapid variations in coastal 
morphology, pollution transport, and hazardous swimming conditions (Huntley and Short 
1992; Ranasinghe et al. 2001). 

Early rip current research stemmed from a need to differentiate them from the 
“undertow” and “riptides” that plagued swimmers and exhausted lifeguards at Venice 
Beach, CA (Shepard 1936). Efforts were made to study their nature and morphological 
importance (Shepard, Emery and Ea Eond 1941). Initial studies were either qualitative or 
crude in their experimental setup, often-employing dyes, surface floats (drift bottles), and 
hours of visual observation (Shepard and Inman 1950). Early investigators concluded 
that a relationship existed between the spacing of the rips and the incident angle and 
groupiness of the incoming wave trains. 

Several theories and empirical equations have been proposed to account for the 
regularity and persistence of rip current spacing. All the theories are based on alongshore 
and cross-shore spatial and temporal variations of radiation stress. The setup of the 
waves in the surfzone does not completely balance the incoming radiation stress 
generated by the waves; therefore, currents are generated in the surfzone to return this 
water offshore (Eonguet-Higgins and Stewart 1964). In nearshore circulation cells, such 
as rip systems, this water will flow downward from regions of high waves (high setup) as 
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feeder currents into the low wave regions where it is returned offshore as a rip current 
(Dalrymple 1975). 

The rip current theories vary by the cause of the alongshore variations in the 
radiation stress. Bowen and Inman (1969) showed that in theory and the laboratory for a 
planar beach that that there is a coupling between an incident wave train and a standing 
edge wave that results in a longshore variation in wave height and radiation stress. This 
in turn results in preferential spacing for the rip currents at alternating antinodes of the 
sub-harmonic standing edge waves of mode n=0. However, this rather special scenario is 
seldom observed in nature. 

Bowen (1969) considered the case of a normally incident uniform wave train on a 
beach with undulating bottom topography. Alongshore variation of radiation stresses are 
created by waves refracting over the undulating bathymetry, from which Bowen obtained 
an analytical solution by linearizing the equations of motion. Circulation cells were 
alculated in which shoreward and offshore flows were of equal streamline spacing. 

Sonu (1972) found that rip current cells were observed with O(60m) spacing at a 
Florida coast for a smooth beach with normally incident waves. The waves tended to 
break in the depressions vice the shoals, in contrast with observations by Shepard (1950) 
off of Scripps Beach. Sonu (1972) considered the cases for alongshore variations in 
radiation stresses due to alongshore variations in wave height offshore and uniform 
waves incident upon an undulating bathymetry. Either leads to the generation of 
nearshore circulation cells. 

Alongshore varying wave momentum flux was generated by Dalrymple (1975) 
using intersecting wave trains. Dalrymple (1975) showed in theory and laboratory, that 
intersecting waves of opposite direction but of the same period generated alongshore- 
varying radiation stresses and a stationary rip system. Fowler and Dalrymple (1990) 
generalized Dalrymple (1975) by considering two groups of directionally narrow random 
waves of opposite but different directions, which resulted in a slowly propagating rip 
current system. 
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Most recently, Reniers et al. (2003b) generalized the work by Fowler and 
Dalrymple (1990) by considering a directionally broad energy spectrum at the offshore 
boundary applied to the non-linear shallow water equations. They show that the groups 
of incident waves arrive at the shoreline with alongshore length scales comparable to the 
rip spacing. The spatially variable energy generates variable radiation stresses that create 
low frequency (~15min period) vortices with alongshore length scales that act to perturb 
the initially alongshore-uniform bathymetry. A positive feedback of the forced vortices 
with the perturbed bathymetry results in an alongshore bar incised by rip channels. The 
spacing is related to the directional spreading and energy of the incoming wave train 
(Figure 1). 

B. BEACH CUSPS 

Rarely is a shoreline longshore uniform. More commonly, due to the complex 
interaction of the waves, currents, and local topography, a coastline is three-dimensional 
and exhibits a variety of features such as offshore bars, rip current cells, and beach cusps. 
Cuspate features ranging from O(lOm) meter spacing on dissipative beaches with low 
wave energy to 0(10,000) meter spacing in the case of the North Carolinian capes are 
prevalent on many beaches (Komar 1998). This study focuses on what are known in the 
literature as “rhythmic topography,” “sand waves,” or “giant cusps,” (1998), (as opposed 
to beach cusps, O(lOm)), of 0(100) spacing that characterize the Southern Monterey Bay 
shoreline and will be referred to simply as cusps for this study (Figure 2). 

The dynamics between beach cusps of various sizes and rip currents is not well 
understood, though they have been documented since the inception of rip current research 
(Shepard, Emery and La Fond 1941). Rip currents were found to exist in the 
embayments of both small and larger cuspate formations of Scripps Beach. In contrast, 
observations by Komar (1971) in the field and the laboratory indicated that the rips 
emerged from the shoals. This is consistent with the relationship between breaking wave 
height and depth found by Mei and Liu (1976), i.e. rhythmic bathymetric contours in the 
nearshore would force the rip from the embayments into the shoals. However, the 
relatively small width of the rip channels compared to the shoals negates the radiation 
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stress imbalance caused by the later wave breaking in the embayments. Consequently, 
rip channels form within the embayments, not the shoals. 

Southern Monterey Bay affords a unique environment for the study of rip currents 
owing to their prevalence created by near normally incident waves impinging on a sandy 
shoreline. Sheltered by headlands to the north and the south, the predominant Pacific 
NW swell from the Aleutians enters the bay as a narrow-banded wave train (Munk and 
Traylor 1947; MacMahan et al. 2002b). The strong refraction over the Monterey Bay 
Canyon tends to bend the wave trains into a near normal approach to the shoreline. This 
refraction also leads to dramatic gradients in wave energy from the sheltered barred 
beach of Del Monte to the energetic, rip-dominated shore of Marina State Park (Figure 
3). It is hypothesized that the spacing of the rip current channels is related to the 
directional spreading of the incoming wave trains. 

Given the difficulties of making in situ nearshore bathymetric measurements, it 
will first be shown that for a cuspate shoreline in the Southern Monterey Bay, there is a 
rip in the center of each cusp. Evidence for the relationship between rip currents and 
embayments in the Southern Monterey Bay is demonstrated in the aerial photograph from 
1984 (Figure 4). Taken on an energetic day at low tide with the sun directly overhead, 
the photo clearly shows rips to be located in the center of each cuspidal embayment. 
Therefore, the more easily measured cusp spacing will be used as a surrogate for rip 
spacing. Beach cusp formation favors normally incident wave trains (Longuet-Higgins 
and Parkins 1962). Given the prevalence of normally incident waves along Southern 
Monterey Bay, cusps and embayments are rhythmic features common to this study area. 
The image only covers a small portion of the coastline, but is representative of the 
Southern Monterey Bay from Sand City until the Salinas River. In the next section, the 
determination of rip current morphology by in situ surveys and remote sensing using 
LIDAR and video imaging (Argus) is described. The wave climate is then described, 
followed by how the rip current spacing is related to the wave climate. 
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II. RIP CURRENT MORPHOLOGY 


A, SURVEYS 

A study of coastal morphology must be able to account for rapid changes in 
bathymetry and beaeh morphology over a range of wave conditions, and especially 
during storm events. Ground truth of an area is necessary to ealibrate the remotely 
sensing LIDAR and Argus systems, as well as any littoral eirculation models. This 
requires periodie personnel-intensive bathymetrie and shoreline surveys. Bathymetries 
and beaeh eontours were obtained using three surveying methods; ATV, personal 
watercraft (PWC), and walking. An ATV rigged with an Ashtee real-time kinematic- 
differential global positioning system (KGPS) receiver surveyed a eross-shore distance 
extending from the maximum run-up at low tide to the dune face, and an alongshore area 
of O(lkm) at each beach survey. To ensure accuracy in the ATV’s position and 
elevation, the Ashtee KGPS was synehronized with the loeal base station at Sand City 
(List and Farris 1999). The KGPS system error when mounted on the ATV is assumed to 
be that of a stationary system; lem in the horizontal and 2om in the vertieal. The same 
KGPS receiver was also mounted on a baekpaek worn by a surveyor in the low tide 
surfzone to a depth of Im and an elevation slightly higher than the closest ATV traek to 
the shore. No error testing has been eondueted yet for the walking survey, but the error is 
estimated to be around lOem in the vertical and horizontal, with variations depending on 
the individual walker (MacMahan Personal Communication 2003). At high tide of the 
same day as the other surveys, a PWC equipped with a KGPS reeeiver and eehosounder 
eovered the area from several hundred meters offshore to about Im of depth. When 
mounted on the PWC, the KGPS error is 6om in both the horizontal and the vertieal. The 
PWC eehosounder has a horizontal spatial resolution of 0.8m and an aeeuraey of 3om 
(MaeMahan 2001). In this manner, each of the surveys overlapped the others. 

Surveys were eondueted over the eourse of the winter season, November to 
February, whieh represents the most energetie eonditions in the Bay, as well as the most 
interesting beaeh morphology. Condueted at around three week intervals during the 
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study, Del Monte, Sand City, and Marina beaches were all surveyed within days of each 
other to make temporal comparisons possible. 

A survey was conducted 16 January 2003 during a time of no wind and a 
moderate swell of H^=\.lm. Low wind and wave energy conditions are ideal to 

conduct the beach surveys. Too large of a significant wave height hinders the PWC’s 
intrusion into the breaking zone and the walker’s cross-shore wading distance. The 
survey of Del Monte Beach shows a low-tide terrace with no rip cells (Figure 5). The 
survey tracks of the ATV, walker, and PWC are traced in the figure, demonstrating the 
overlap between techniques. 

The survey of Sand City Beach on the same day (Figure 6) shows a low tide 
terrace incised by rip channels with spacing of 0(150m), similar to results during the 
RIPEX experiment at this same beach (MacMahan 2003a). The GPS tracks indicate that 
full coverage of the survey area was not achieved. The lower-most shoal in Figure 6 is 
generated by the Delaunay triangulation that drives MATLAB’s gridding function 
(Barber et al. 1993), but the lack of sampling in that region causes those contours to be 
suspect. It is perhaps easier to distinguish the bathymetric features with the contour 
image (Figure 7); however, the cuspate features of the beach are more easily discernible 
in the colored plots (Figure 6). The cusps are of the same order (150m) and occur 
concomitantly with the rip channels. This supports the findings of Bowen and Inman 
(1969), Komar and McDougal (1988), Revell et al. (2002), and Sallenger et al. (2002) 
that the rip current channels occur in the embayments of cuspidal shorelines. 

The Marina State Beach was surveyed on 16 January 2003 (Figure 8) and shows a 
well-resolved shoal at y = -100m and a cuspidal shoreline. Qualitatively, the cuspidal 
spacing is O(350m) with the embayments at the mouths of the rip channels. The beach 
and inter-tidal area was well sampled to a depth of Im, but the sparse PWC echosounder 
data makes the analysis of the rip spacing uncertain. Thus, the rip channel spacing is 
determined from the location of the rip channels. 

Initial surveys suggest that to ensure full coverage of the beach, it would be 
beneficial to the surveyors if visual markers were erected. Often the PWC, ATV, and 
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walker did not cover the same alongshore areas; this led to gross interpolations by the 
Delauney bi-linear triangulation method at the upper and lower alongshore bounds 
(Figure 7). The surveys suggest that a hatched diagonal pattern might be adopted by the 
PWC and a few alongshore transects be added to the walker’s survey. This would better 
resolve some of the embayments and beach cusps, as well as avoid many artifacts in the 
estimated topography. Also, unless the longshore bar is being sampled, the surfzone 
features of interest generally have alongshore variations, which are difficult to resolve 
with shore normal PWC transects. For example, offshore transverse bars at Duck, NC 
were not observed until captured by time-exposed video imagery because previous 
surveys of cross-shore transects had poor alongshore resolution (Konicki and Holman 
2000 ). 

B, LIDAR 

As mentioned in the introduction, it was qualitatively found that a rip current is 
concomitant with each cusp embayment for the Southern Monterey Bay (Figure 4). 
Given this relationship, we consider cusp spacing, , to be a proxy for rip spacing, . 

Beach surveys are used here to determine cusp spacing. 

Light Detection and Ranging (LIDAR) surveys were conducted before and after 
the 1997-1998 El Nino winter to measure erosion along the southern coast of the 
Monterey Bay (Thornton et al. 2003). Determining beach cusp dimensions along the 
southern coast of the Monterey Bay is a byproduct of that 1997-98 study. The LIDAR 
imaging was part of a larger effort by the United States Geological Survey (USGS), 
NASA, and NOAA to map the erosional response of the western coast of the United 
States during that El Nino (Sallenger et al. 2002). The EIDAR survey was conducted 
using NASA’s Airborne Thematic Mapper (ATM) mounted underneath a plane (Eigure 
9) and calibrated for 3-D positional accuracy as described in Brock et al. (2002). The 
ATM has a swath width of 350m given an aircraft altitude of 700m, a horizontal point 
resolution of L4m (Brock et al. 2002), and a vertical RMS error of less than 15cm 
(Sallenger et al. 2003). 
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Post-processed LIDAR images generated by the NASA’s LaserMap program are 
loaded as geotiff files into the AreView Geographie Information System (GIS) suite 
(Broek et al. 2002). Images are overlain on Digital Orthophoto Quadrangles with a final 
output displayed in UTM eoordinates with NAD-83 as the horizontal datum (Figure 10). 
The underlying images were taken 21 August 1998, after the LIDAR survey. A shoreline 
is defined by MLLW, with watered areas excluded from the output by use of an edit 
mask. Though over-sampling the data, a grid spacing of 1 m was chosen to represent the 
variable ATM sampling rate (Thornton et al. 2003). It was deemed desirable to resolve 
beaeh eusps O(30m) in the analyses, whieh are diffieult to resolve both near the waterline 
(0 m contour) and near the toe of the dune (5m eontour). Therefore, the two-meter 
eontour was ehosen to represent the shoreline. 

The 2m eontour, plotted in UTM coordinates, was rotated 60° to provide a 
monotonie X-eoordinate. A 2”‘*-order polynomial was fitted, in a least square sense, to 
the data set and subtracted to remove the curvature of the Southern Monterey Bay (Figure 
11). Some distortion is present in the vieinity of Del Monte Beaeh and the Salinas River, 
identifiable by the positive trend in y-axis deviation from zero, moving left and right 
from the eenter of the x-axis, respeetively. The resulting error in shoreline length is less 
than 1.5 pereent in these regions. To remove the distortion, the data were then high-pass 
filtered, removing variations greater than 1000m scale. The resulting shoreline contains 
cuspate features oscillating about a straight shore (Figure 12; Figure 13). Qualitatively, 
the spaeing of the cusp features increases in length and in width from south to north after 
the El Nino winter. 

Wavenumber speetra (m ') were generated for both surveys (Figure 14; Figure 
15) to examine the lengths of the euspate features. The speetra are truneated for 
wavenumbers greater than 0.05 (m') as energy at these wavenumbers are features too 
small to be cuspate. The speetra were caleulated by subdividing the full reeord into 
1200m lengths resulting in twenty-two degrees of freedom for both the 1997 and 1998 
surveys. The signals from both years indieate a broad peak with a peak eorresponding to 
a 340m eusp width. The energy speetra show that for both years the dominant 

oscillations are similar to the spacing of rip currents. The broad spectral peak about the 
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value of 0.0029m ' (^c«^p~340m) is due to speetral smearing owning to the 
significant variability in cusp spacing. Beach cusps are evident in the 1998 LIDAR study 
as indicated by the broad spectral bump centered on wavenumber 0.021m ' 
~40m). 

To resolve the variable cusp spacing alongshore, the data were segmented into 
multiple bins approximately 1200m in length. Each segment contains complete cusp 
wavelengths to reduce leakage in the spectral analysis. Segmentation decreases the 
degrees of freedom available from the length of record, but the 1200m record length was 
chosen to measure variability and give a reasonable cusp statistic (5-10 cusps per 
segment). 

The average cusp spacing for each bin, , was calculated as the peak of the 

wavenumber spectrum (Figure 16). Prior to the winter, the cusp spacing initially 
increases and then is more or less constant. After the El Nino winter, there is an increase 
in cusp lengths proceeding north toward the Salinas River. The increase in cusp spacing 
moving up the coast, and hence the increase in rip spacing, is consistent with both our 
observations and hypothesis. 

The average horizontal cusp excursion, 

where H is the horizontal excursion measured from the center of the embayment to a 
point parallel to the horn for each segment (Figure 17). Two trends are noticeable from a 
comparison of the cusp excursions of the two surveys. First, the excursions grow in size 
for both years as one proceeds from Monterey northward along the beach toward a more 
energetic wave climate. Second, the cusp excursions are noticeably greater in magnitude 
for the 1998 survey. 

A steepness parameter is calculated as a measure of the horizontal excursion of 
the cusp over the alongshore cusp length; 
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(2.1.2) 


m^=riJX^ 

The steepness is fairly eonstant for fall 1997, until alongshore distanee 6000m, whereas 
the spring 1998 survey has much greater cusp steepness values along the entire length of 
the surveyed shoreline (Figure 18). This probably reflects the increase in rf^ due to the 
scour and increased erosion deepening the cusp embayments during the El Nino storms. 

Another LIDAR survey by the United States Geological Survey is planned for the 
Monterey Bay during the summer of 2003. However, the survey will only provide a 
snapshot of the rip spacing in the Bay. Besides cost and time constraints, LIDAR does 
not provide insight into the temporal variability of morphology. A need exists for auto- 
archived and ongoing remotely sensed time-series of nearshore processes. The Argus 
camera system fills that gap. 

C. ARGUS 

The Argus video imaging system is a means to obtain long term and beach 
information without the cost of traditional in-situ instrumentation. The Argus system 
does not survey to the same accuracy as conventional techniques, e.g. KGPS walking 
transects, but it is attractive for its cost, immunity to weathering, and temporal sampling 
ability (Stockdon and Holman 2000). 

So far, only video images at Marina State Beach are processed on a regular basis. 
The Marina Argus station is an array of five video cameras mounted approximately 15m 
above the shoreline. This area of coverage is greater in the alongshore than that of the 
GPS surveys to date at Marina. The Sony 640x480 pixel cameras have an image 
resolution dependent on the distance to the object of interest, but typically are of the same 
order as GPS survey resolutions (1997). They are housed within standard, 
environmentally resistant surveillance camera housings. The protective glass window 
has a tendency to fog up or get dirty at this location due to its proximity to the beach and 
its location at the edge of the dune. The Marina cameras were initially lined up in late 
October, but were vandalized by bending them out of position and smearing them with 
surf-wax between 02NOV02 and 06NOV02, 1900GMT. They were rectified and 
working properly until 09FEB03, 1900GMT, at which point they were again vandalized 
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and knocked out of alignment. They were realigned on 20 February 2003. The five 
eameras are angled to eover more than 180 degrees of nearshore area with minimal 
overlap, resulting in a colleetive Field of View of 2 km alongshore and 1/2 km aeross 
shore. 

Eaeh eamera generates a snapshot of the beaeh on the hour on a daily eyele, 
resulting in images from after dark being stored in the database as well. Also, a ten- 
minute time exposure (timex) is created eaeh hour to show the general trends of the surf 
at that time. This exposure length was found to be ideal for Atlantie wave eonditions and 
has been adopted in the Monterey Bay (Holman and Lippmann 1987). Finally, highly 
dynamic areas of the surf zone are shown in a varianee image derived from the timex. 
43200 stored images are generated over the course of the year. Needless to say, one must 
be highly seleetive of whieh images will be analyzed, or even viewed. Many may be 
eulled based on the times listed, e.g. those at night, but others must be loaded and then 
visually eheeked for unsuitably due to foggy/low visibility eonditions, salt or grime on 
the lens, sun-indueed glare, misalignment of the eamera, ete. 

Ground Control Points (GCPs) are required to solve the equations needed to 
reetify an Argus image so that an objeet is known absolutely in both image and world 
eoordinates (Holland et al. 1997). Preferably, GCPs of opportunity are used to save on 
eost and time. They might inelude a pier piling, tall building, fenee post, or any other 
objeet that will remain in a relatively fixed loeation during the course of the experiment. 
The Marina site has a fenee immediately in front of the eameras that provides several 
GCPs; however, its elose proximity to the array means that the aeeuraey of its location 
must be great, or severe distortion will occur at the far edges of the eamera’s viewing 
field. There are a few other GPC’s, including the horizon, that ean be used in eaeh of the 
five images generated by the array. Despite these natural validation points, a few 
manufaetured GCPs are brought in during reetifieation to ensure the proeess is aeeurate 
and points are seattered throughout the images. Self-employed GCPs are typieally large 
white spheres that are placed in the eameras’ viewing window, have their GPS 
eoordinates reeorded, and are identified in a grayseale image by their large eontrast with 
the beach and ocean (Holland et al. 1997). Rectification allows the scene captured by a 
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camera’s lens to be transformed from image coordinates (U and V) into more familiar 
world eoordinates (X and Y) (Konieki and Holman 2000). It standardizes the different 
views presented by each camera in the array, and allows their eollective image intensity 
maps to be smoothly merged into one plan view of the beach (Figure 19). For a three 
eamera system at elevation 10m and near to the survey area, Holland et al. (1997) 
obtained world coordinate errors of 0(lm) for the rectified images. This order is the 
same as the error by traditional surveys (Salinger 2003) and is extended to the five 
camera Argus system at Marina. 

Previous experiments have successfully compared traditional in situ 
morphological surveys with Argus timex imagery. Konieki and Holman (2000) 
confirmed the location of transverse sand bars at Duck, NC using both the Coastal 
Research Amphibious Buggy (CRAB) and an Argus array. However, accurately 
identifying individual rip occurrences from these images is still a naseent technique. 

Ranasinghe et al. (1999) used zero upcross/downcross analysis to locate rips at 
Palm Beach, Australia. Those rips were easily distinguishable by visual observation at 
the shoreline as dark patches of water with little foam or wave breaking, and 
consequently as troughs in the longshore transects of pixel intensity generated from the 
Argus time exposures. This teehnique is not adequate for Marina State Beach, as the rips 
tend to be obscured by resident foam and surfaetants on the water’s surfaee and therefore 
present smaller longshore pixel intensity gradients in the images (Symonds Personal 
Communication 2003). Also, the Marina camera is only 15m above the survey area, 
while the Palm Argus setup commands a view from approximately 100m up on a bluff 
overlooking the survey area. A view from a higher incidence angle would allow for a 
greater field of view and capture more rips in the time exposures. 

Another method employed by nearshore researchers is to study the lagged 

autocorrelations function (ACF) of the alongshore intensity transects. Quantitatively, the 

lag of the local maximum of the eorrelation coefficient is noted as the rip spacing for that 

imaged area (Ranasinghe et al. 1999). Peaks in that function at a 95% confidence 

interval were looked for by Ranasinghe et al. (1999) in their Palm Beach, Australia 

analysis. Unfortunately, the rips laeked periodieity in their spacing and the ACF 
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generated no statistically significant peaks. Given the gradient in rip spacing 
demonstrated by the morphology and LIDAR surveys, we would be ill advised to use this 
method in the Monterey Bay. 

Though not as rigorous, this study estimates from visual analysis of a large 

sample of timex images at Marina during the months of the GPS beach morphology 
surveys. Preferential wave breaking over a shoal is represented by a large pixel gradient 
in that area of the timex image. Gradients in bathymetry are inferred from these pixel 
gradients. Rip currents are clearly visible in the rectified Argus images of Marina Beach 
as areas of low intensity pixels extending seaward with surrounding high intensity pixels 
of shoal-breaking waves. Close attention was paid to 16 January 2003 and 14 February 
2003, dates when all three beach survey methods worked in unison. Unfortunately, the 
cameras were jarred out of alignment just prior to 14 February 2003 and were not 
geometrically realigned by that time. Consequently, images from 8 February 2003 and 
14 February 2003 were used as proxies, with consequent reservations about temporal and 
spatial applicability (MacMahan 2003a) for the actual survey date. Average spacing 

during the study was of O(300m) at Marina Beach. The method employed to measure rip 
spacing via video imagery at Marina is accomplished by hand analysis of images 
generated once weekly from November 2002 to April 2003. This qualitative assessment 
of rip spacing over the energetic winter season gives estimates of O(300m). This is 

in agreement with our morphological surveys of the area (Figure 8), as well as the 
energy-plots of that bin from the LIDAR survey (Figures 14 and 15). 

Contoured images from the beach surveys are overlaid on their respective Argus 
images (Figure 20). Qualitatively, a strong correlation is present between the breaker¬ 
line in the Argus image (represented by the line of white pixels just offshore), and the 
narrowing of the contours in the GPS survey images (rhythmic longshore bars). The 
cusps identified by the walking survey at Marina on 16JAN2003 are almost perfectly 
overlaid with the white pixel undulations between 50-100m across shore. The rhythmic 
pixel pattern in the Argus image repeats at 0(350m). Ground truthing of the shore 
normal alternating streaks of white and dark pixels as rip channels and connected bars is 
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difficult given the incomplete coverage of the PWC survey during any of the Marina 
surveys. Though not eoneurrent, the Argus timex of 08FEB2003 and the PWC transeet 
of 14FEB2003 both indieate the presence of a strong rip eurrent at approximately x = 0m 
alongshore (Figure 20; Figure 8). Unfortunately, the PWC surveys only eovered 20% of 
the Argus field of view, thereby preventing any other features from being validated. No 
walking or ATV survey was arehived for that date, so a verification of Argus-imaged 
beach cusps is not available. 
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III. WAVE CLIMATE 


A, WAVE ENERGY AND DIRECTION 

To test the hypothesis that the rip current spacing is a function of wave energy 
and or wave directional spread, directional wave parameters are compared between 
measurements obtained from two Datawell Directional-Waverider buoys at Marina State 
Beach and Sand City Beach, both in 18m of depth (Figure 21). These two sites are 
separated by ~6km. To compare the variation in wave parameters along the shoreline of 
the Southern Monterey Bay, a time period of yeardays 120 to 188, 2002 was chosen for 
analysis when both waverider buoys were operating. These measurements are used to 
test both a refraction and nearshore circulation model. 

The energy of an incoming wave train is described by the directional spectrum; 

s{f,e) = s{f)D{f,e) (3.1.1) 

where 0 signifies the arbitrary incidence angle of a wave train. The directional spectrum 
can be represented as a Fourier series in direction at a particular frequency; 

Dif, 0) = - I + X cos(n^) + sin(n6')} 

«=i J j 2) 

= -^ + aj cos ^ + hj sin ^ + aj cos 20 + sin 20 

Longuet-Higgins et al. (1963) demonstrated that the first four Fourier coefficients could 
be derived from the spectra and cross-spectra of the pitch, roll, and heave of a buoy or 
equivalently, the pressure and horizontal velocity components (p,u,v). 

<^i(/) = J d0cos0S(0;f) 

c,Af) 
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(3.1.4) 


b,{f)=r ddsm0S{d;f) 
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(3.1.5) 


(3.1.6) 


Both the buoy and ADCP data are arehived with energy and these low-order eoeffieients 
at each frequency bin. 

The waves at the two sites are compared using parameters for height, direction, 
and energy flux. Significant wave height is calculated: 

I r 0.25Hz 

where the frequency integral is taken across the wind-swell energy band (0.04-0.25 Hz), 
is similar at the two beaches for this springtime series, with Marina having slightly 

higher waves on average (H^ =1.3m) compared to Sand City,//^ =l.lm (Figure 22). 

The exceptions are at yeardays 142 and 152. In those instances, the significant wave 
height is much greater at Marina. Yearday 160 shows an event in which the wave energy 
is greater at Sand City. This seems to be an anomaly for this time series and is not 
repeated elsewhere. 

The peak period, , is calculated from the energy spectra as the inverse of the 
frequency containing the maximum amount of energy: 

(3-C8) 
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The peak periods , as are expected, do not show much variation between the two 
locations (Figure 23). 


The mean direction of the peak energy is defined in terms of (3.1.3) and (3.1.4): 


0 peak — tan 




(3.1.9) 


where the peak represents the frequency of maximum incident energy at the sensor, 
which is valid only if S{f ) is narrow-banded (Herbers et al. 1999), as is most common 
in the Monterey Bay (MacMahan et al. 2003). Refraction within the bay tends to result 
in waves approaching the shoreline shore-normal. Marina has 6peak = T and Sand City 

Opeak = 2.5°, indicating that the wave climate at Marina is not as strongly refracted as that 
at Sand City (Figure 24). 

Peak and bulk directional spreading parameters are calculated from the directional 
coefficients: 


<ypeak=m-alpeak^hl^eakTT O-LW) 

( 3 . 1 . 11 ) 


which are also strictly valid only for narrow-banded S{f ,6) . The spreading coefficients 

represent the half-width location of S{6,f ) when plotted against Opeak and Obuik, and 
indicate the distribution of energy about the mean direction (Figure 25). 

The measured peak directional spreading parameter (3.1.10) was the same for 
both Sand City and Marina, =14°. The bulk directional spreading parameter was 

slightly larger at Marina, = 20°, compared with =18° at Sand City. The larger 
a^uik Marina suggests greater directional bandwidth at higher frequencies, as this site is 
more open to waves and wind from offshore. 

Alongshore sediment transport is a function of the alongshore energy flux, . 

Hence, energy flux is used to parameterize morphology, which is a manifestation of 
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gradients in sediment transport. Both along- and eross-shore energy fluxes are 
caleulated. The measured energy flux is obtained from the direetional speetrum: 

F^ = pg\'^dfc^{f)\yecoses{f,e) ( 3 . 1 . 12 ) 

where p is the density of seawater, g is the aeeeleration of gravity, / and represent 
the wave frequeneies to be integrated aeross, and c^is the group veloeity. Using(3.1.3) 
and (3.1.4), the eross-shore bulk energy flux is: 

r^.lSHz 

dfc^{f)a, (3.1.13) 

and the alongshore bulk energy flux is: 

r^.lSHz 

(3.1.14) 

The signs of both the flux parameters indieate the direetion of bulk transport of 
the wave energy at each site. Negative values indicate onshore transport, while 

positive F^, values represent an energy flux directed northwards toward the Salinas 
River. 

Both Sand City and Marina (Figure 26) exhibit strong onshore transport of 
energy. There is a noticeable gradient between the two sites with Marina having an 

average cross-shore flux of Fx = 33000(J/ m^) and Sand City Fx = 28000(J/ m^). This 
gradient is not surprising given the significant wave height differences between the two 
sites. 

There is significantly more variation in the direction of the longshore energy 
fluxes for both Sand City and Marina (Figure 27). Marina measured Fy = 3250(J/ nd), 

while Sand City measured Fy = 440(J/ nd) . These values are a result of the increased 
incidence angle at Marina compared with Sand City. 
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B, COMPARISON WITH REFRACTION MODEL 

Directional wave spectra are calculated by refracting the directional wave spectra 
measured at NOAA buoy 46042 (Figure 28) located in deep water outside Monterey Bay 
(O’Reilly and Guza, 1998). From the refracted directional spectra, time series of and 

Opeak are generated for locations corresponding to the two Datawell buoys. The model 
bulk energy is calculated for the same swell-wind frequency bands as the in-situ buoys 
use. 

The measured and refracted model values are compared at Sand City and 

Marina (Figures 29 and 30). Linear regressions of the modeled and in-situ data show the 
poor correlation at both of the locations (Figures 31 and 32). The refraction model 
underestimates the significant wave heights by about 60% at Marina and 50% at Sand 
City. 
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IV. RIP CURRENT SPACING RELATED TO WAVE CLIMATE 


A. LIDAR 

No direct comparison of the LIDAR 2m contour analysis and wave climate is 
possible as no working directional wave buoy was deployed within the Monterey Bay 
during that time, and NOAA’s offshore directional wave buoy just outside of the Bay was 
inoperable that winter. The LIDAR imagery only contains data from pre- and post-El 
Nino events; therefore, it might be premature to extend these findings to non-El Nino 
winter shoreline morphology in the Monterey Bay (Sallenger et al. 2002). 

EIDAR was also applied to the analysis of the 1997-98 El Nino-induced erosion 
in the Netarts Eittoral Cell, Oregon (Revell, Komar and Sallenger 2002). They found that 
greatest amount of erosion in the dunes of the Netarts Cell occurred in the rip current 
embayments. Extension of this finding to the Monterey Bay would support our 
hypothesis that the rip currents are located at the center of the embayments. 

The wavelength of cusps increased to the north, so they are not spatially 
homogenous. Using another analysis technique, such as wavelets, instead of Eourier 
transforms might increase the degrees of freedom for our shortened segments and hence 
the accuracy of our results. That method should be considered for future research. 

B. ARGUS 

Argus allows the determination of morphological changes in response to waves 
on an hourly basis. A current proposal by the Naval Postgraduate School Oceanography 
Department calls for the replacement of the existing camera housings with improved 
versions equipped with automated windshield wipers at Marina. Additionally, cameras 
will be added at Sand City and the MISO site to allow for better quantification of the 
morphological variation along the shoreline of the Southern Monterey Bay over time. 
Combined with the improved GPS surveys and contouring methods (objective analysis) 
at the MISO site. Sand City, and Marina scheduled for this summer, a more quantitative 
analysis of the rip spacing from the Argus imagery will be possible. 
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C. WAVE CLIMATE 

Hopefully, future analysis will allow wave elimate and morphologieal data to be 
eompared for the same time period. It is presumptuous to draw eonelusive results and 
parameter relationships for rip spaeing from wave elimates of a different season. As an 
example, for the Marina site, elimatologieal wave height measurements over 5.5 years in 
the early 1990’s show an average differenee of one meter in signifieant wave height 
between the energetie winter season and the more ealm summer season. Our 
morphologieal beaeh surveys ehronieled the monthly ehanges in bathymetry at the three 
sites during the winter of 2002-03, whieh was also an El Nino winter. The LIDAR 
surveys only represented snapshots of the eonditions at two short intervals in Oetober 
1997 and April 1998, respeetively. No details were doeumented during the 1997-98 El 
Nino winter. Wave statisties were analyzed for yeardays 122-188 of 2002 beeause this 
was when both of the Naval Postgraduate Sehools direetional wave-riding buoys were 
operational. 

D, REFRACTION MODEL 

Signifieant efforts have been undertaken to model the direetional distribution of 
nearshore wave trains (Munk and Traylor 1947; Abreu et al. 1992; O’Reilly and Guza 
1991,1998). Data from wave buoys is erueial for in-situ information and ground truthing 
of modeled data; however, many eoastal areas, espeeially those alien to our shores but of 
interest to the Navy and the Marine Corps, do not have suffieient buoy eoverage to allow 
for preeise interpolation of loeal wave elimates and nearshore eireulation. 

O’Reilly’s pure refraetion wave model was initialized and run for the same time 
period as the wave buoy measurements. The model under-predieted the wave heights 
eompared with the buoy measurements by 0(50%) at Sand City and Marina. Efforts are 
being made to improve the refraetion model. The hope is that onee the refraetion model 
is verified, it ean be used to prediet the wave heights, energy flux, and direetional 
bandwidth along the entire shoreline of the Southern Monterey Bay. In addition, it is 
hoped that the refraetion model ean extrapolate NOAA buoy measurements of wave 
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climate off the Farallon Islands (the nearest wave buoy to the north) to Monterey Bay 
during the time of the 1997-98 LIDAR measurements. 

E. MORPHODYNAMIC MODEL 

Direetional spreading is a parameter that has only reeently been ineluded in 
nearshore eireulation models. Renters et al. (2003) ran ealeulations using the 
morphodynamie model Delft 3-D for various values of direetional spreading with mean 
angle normally ineident and holding wave height eonstant, = \m (Figure 1). The 

absenee of direetional spreading, 0 -^= 0 °, led to a minimum spaeing of the rip eurrents, 
eonsistent with a self-organizing generation meehanism. A small amount of <jg 
transformed the beaeh into a quasi-foreed system with preferred rip spaeing governed by 
the ineoming wave groups (Figure 1). Further inereasing of cjg eaused the rip spaeing to 
deerease as the alongshore length seales of the wave groups grew smaller (Figure 33). 

Given that both Marina, (Jbuik=2Q°, and Sand City, (Jbuik=\T, fall on the 
shallow negative slope of Renters’ model (Figure 33), there should not be a large 
differenee in their rip spaeing for the same wave heights. The nearshore model prediets a 
spaeing of approximately 0(180m) for the two sites for H^=lm . These spaeings are 

signifieantly less than the aetual measured for Marina, = 340m, and slightly less 

for Sand City, =150m. Renters (Personal Communieation 2003) has run the Delft 

3-D model and found that the rip spaeing inereases with inereased wave height, 
eonsistent with the observations in Southern Monterey Bay. Future use of the model 
should foeus on varying the ineident energy of the wave fields. 
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V. CONCLUSIONS 


Rip current spacing in the Southern Monterey Bay is compared with the 
alongshore gradient in wave climate to determine if a relationship existed. It is 
hypothesized that the energy and directional spreading of the wave trains determined the 
spacing of the bay’s rip ehannels. Monthly beach morphology surveys during 2002, 
LIDAR surveys from the 1997-98 El Nino winter, and Argus time exposures from the 
winter and spring of 2002 are used to quantify the width of the rips both spatially along 
the eoast and temporally throughout these seasons. Significant wave height and 
directional spreading parameters are eompared using measurements obtained for 
direetional wave-rider buoys at Sand City and Marina for yeardays 120 to 188 for 2002. 
The hypothesis is tested against the in-situ data, and then against the O’Reilly Refraction 
model and the Delft 3-D morphology model. Direetional spreading does not vary 
substantially between Sand City and Marina for the climatology data, and therefore does 
not appear to be signifieant for determining rip eurrent spacing in Southern Monterey 
Bay. Signifieant wave height and energy flux should be researehed further to assess their 
relation to the rips. 

The sites of Sand City and Marina State Beaeh along the Southern Monterey Bay 
eoastline withD 6 km separation between them (Figure 21) allow for general eonclusions 
regarding the gradient of wave parameters and rip current spaeing along the southern 
eoast of the Bay. Rip spacing does indeed inerease along the Southern Monterey Bay 
towards the Salinas River. This hypothesis is validated with qualitative rip spacing 
measurements from beaeh morphology and Argus surveys, as well as quantitative 
LIDAR-generated measurements of the euspate oscillations of the Southern Monterey 
Bay as a proxy for rip spacing. A spacing of 0(150m) is found for Sand City, while the 
spacing for Marina is O(340m). 

It is hypothesized that eaeh eusp embayment eontains a rip channel for the 
Southern Monterey Bay. The order sealing of the alongshore wavenumber, 
= 2.9x1 ' = 340m), from the LIDAR 1997-1998 eusp analysis is in 
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agreement with the measured rip widths, = 340m, for the Southern Monterey Bay as 

measured by beaeh morphology surveys and Argus imagery during winter 2002-2003. 
Therefore, eusp spaeing is a good surrogate for rip eurrent spaeing. 

Argus images are eapable of remotely measuring rip eurrent spaeing over large 
spatial and temporal seales. Qualitative measurements of rip eurrent spaeing during 
winter 2002-2003 at Marina are found to be O(350m), in agreement with beaeh 
morphology surveys during the same period at that beaeh. They are obtained by inferring 
bathymetry gradients from pixel gradients in the images. Offshore features at Marina are 
sometimes diffleult to distinguish, either due to dirty eamera lenses and a low eamera 
elevation above the beaeh, or to resident foam and surfaetants smoothing the pixel 
gradient in the time exposures. 

Peak Direetional Spreading, apeak = 14°, does not vary between Sand City and 
Marina during the study, while bulk direetional spreading has a small variation with 
Marina a bulk =20° and Sand City a bulk =18°. Direetional spreading does not appear to 
be the dominant faetor in the rip eurrent spaeing in the Southern Monterey Bay. The 
gradient of both signifieant wave height, and eross-shore wave energy flux between Sand 
City and Marina should be eonsidered in relation to the rip spaeing and nearshore 
morphology for these loeations in future studies and predietive models. 

O’Reilly’s wave refraetion model under-prediets by 0(50%) the signifieant wave 
heights ealeulated by the direetional wave buoys at both Sand City and Marina. 
Therefore, O’Reilly’s model eould not be used to extrapolate the wave elimate along the 
shoreline of Southern Monterey Bay. 

The diffleulty of this study is that the wave climate that generated the cusps and 
rips of varied lengths is not known. The cusps could have been generated by a single 
extreme storm having particular characteristics, and then persisted for a long duration. 
Therefore, we cannot make definitive conclusions regarding our hypothesis. 

It has been demonstrated that the tools to answer the hypothesis are now in place, 
including the Argus system to examine temporal process changes, wave measurements. 
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and periodic in situ surveys for verification. In addition, frequent 2m contour surveys are 
being conducted by driving along the beach with a KGPS mounted on the ATV. 
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FIGURES 


Figure 1. 
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Example of rip currents and underlying bathymetry of a barred beach incised 
by rip channels generated by a research version of the Delft 3-D morphodynamic 
model (Renters et al. 2003). Erosive areas are shown in red while areas of 
accretion are in blue. Waves on the outer boundary are directionally broad with 

= \m and 6 mean = 0°. The initial bathymetry is an alongshore uniform barred 
profde. 
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Shoreline and sediment organization oseillations as a funetion of frequeney 
and alongshore wavenumber . The euspate features of interest vary between 

one day-^one month temporally, and 100m-^500m spatially (from Komar 1998). 
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Yearly Mean Wave Energy m the Southern Monterey Bay 



0 2 4 6 8 10 12 14 16 18 

Distance Alongshore (km) 

Monterey Salinas River 

Figure 3. Mean wave energy derived from twenty years of offshore elimatology 
refracted to the 10m contour, showing a large alongshore gradient in wave climate 
exists in the Southern Monterey Bay due to the shadowing effects of Point Pinos 
at Monterey Wharf Number 2 (0 km) and the wave focusing at the Salinas River 
(18km) from refraction by the Monterey Canyon in the center of the Bay. 


30 






Figure 4. 



Aerial photo of Southern Monterey Bay taken at low tide on April 1984 with 
energetie wave eonditions showing large beaeh eusps assoeiated with rip eurrent 
eells. The eusps and rips are both of 0(100-400m) (Thornton et al. 2003). 
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Figure 5. Morphological survey of Del Monte Beach, Monterey Bay, CA using 
kinematic GPS and an echosounding transducer. The black lines represent, from 
left to right, the tracks of the ATV, walker, and PWC. The morphology is a 
shallow sloped barred beach with no visible rip channels. 
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Figure 6. 


Sand City (Origin 121.8533W, 36.6163E): J011603, 0011703, W011503 
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Morphological survey of Sand City Beach, Monterey Bay, CA using 
kinematic GPS and an echosound transducer. The black lines represent, from left 
to right, the tracks of the ATV, walker, and PWC. A system of shoals and 
embayments is clearly visible with spacing O(150m). 
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Figure 7. 


Sard City (O'igin 121.8533W, 36.6163E): J()11603, 0311703, WOIISOS 
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The survey results from 16 January 2003 at Sand City, Monterey Bay, CA 
using vertical contours. The rip current embayments are clearly visible at 
y = -50m and y = 100m. 
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Marina (Origin 121.8091W, 36.6998E): J011603, G011603, W011603 



Figure 8. A cuspate shoreline with rip current circulation cells of O(300m) is evident at 
Marina State Beach, Monterey Bay, CA from this three-part GPS/echosounder 
survey on 16JAN2003. Rip currents are located O(200m) apart at Y=200m and 
Y=-200m. 
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Figure 9. 
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A typical Airborne Topographical Mapper flight: surveys conducted October 
1997 and April 1998. Horizontal point resolution is 0(1.4m) and vertical RMS 
resolution is 0(15cm). 
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Figure 10. 



Example of a proeessed 2m contour generated by NASA’s Airborne Thematic 
Mapper overlain on a difference survey, with both overlain on an aerial photgraph 
of the Southern Monterey Bay taken near the time of the 1998 LIDAR survey. 
The black line represents the 2m contour in 1997, with 1998 in red. This contour 
was chosen specifically to resolve beach cusps 0(3Om) and larger features. 


37 





meters 


Figure 11. 



Raw signal of the 1997 LIDAR data, plotted in UTM coordinates. The curve 
was rotated 60° to insure monotonicity along the x-axis. 
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Figure 12. 


Figure 13. 
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meters 

2m eontour for Southern Monterey Bay from the 1997 LIDAR survey, after 
subtracting a 2“‘^-order polynomial to remove the curvature of the coastline and 
high-pass filtering to remove oscillations greater than 1000m. 



2m contour for Southern Monterey Bay from the 1998 LIDAR survey, after 
subtracting a 2“‘^-order polynomial to remove the curvature of the coastline and 
high-pass filtering to remove oscillations greater than 1000m. 
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Figure 14. 
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Wavenumber spectrum for the 1997 LIDAR survey of the Southern Monterey 
Bay shoreline. Alongshore wavenumber peaks correspond to cusp spacing 200- 
340m. 
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Figure 15. 
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Wavenumber spectrum for the 1998 LIDAR survey of the Southern Monterey 
Bay shoreline. The peak at K^, =0.003(l/m) corresponds to a cusp spacing of 

= 340m. 
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Figure 16. 



Cusp spacing, averaged at discreet bins O(lkm) from LIDAR surveys of 
Southern Monterey Bay for 1997 and 1998. 
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Figure 17. 



Horizontal Cusp Excursion of the Southern Monterey Bay. The values are 
averaged over bin sizes O(lkm) and are derived from signal analysis of the 1997 
and 1998 LIDAR surveys of the Monterey Bay. 
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Cusp Steepness 


Figure 18. 



Cusp steepness, rjj averaged over O(lkm) bins for the length of the 1997 
and 1998 LIDAR surveys of the Southern Monterey Bay. 


44 








Figure 19. 


ARGUS Timex: 10603 1120 



Along Shore (m) 


An Argus merged and reetified timex from the five eamera array at Marina, 
Monterey Bay, CA from 06 January 2003, showing the cuspate shoreline and 
underlying bathymetry. 
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Figure 20. 


ARGUS Timex: 011603 1200 with GPS Survey: Marina (Origin 121.8091W, 36.6998E): J011603, G011603, W011603 



Along Shore (m) 

An Argus timex from 16 January, 2003 overlaid with eontouring from the 16 
January beaeh morphology surveying of Marina State Beaeh, Monterey Bay, CA. 
A well-defined rip channel is present in both surveys at x = 0m, though the 
alongshore swath width of the PWC limits validation of the cuspate features 
present in the timex. 
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Figure 21. 



Location of the wave climate sensors in the Southern Monterey Bay. 
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Marina 




Yearday 2002 


Figure 22. Significant wave height (m) plotted against yearday in 2002 from the 
directional wave buoys offshore of Marina and Sand City, Monterey Bay, USA. 
Mean values are shown in red. 
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Marina 


Figure 23. 



Sand City 



Yearday 2002 


Peak period (see) of ineident wave trains plotted against yearday in 2002 from 
the direetional wave buoys offshore of Marina and Sand City, Monterey Bay, 
USA. Mean values are shown in red. 
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Figure 24. 
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Mean direetion (deg) of the incident wave trains at peak period plotted against 
yearday in 2002 from the directional wave buoys offshore of Marina and Sand 
City, Monterey Bay, USA. Mean values are shown in red. 
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Marina 


Figure 25. 
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Bulk directional spreading, , plotted against yearday in 2002 from the 

moments of the directional wave buoys offshore of Marina and Sand City, 
Monterey, USA in 18m of depth. Mean values are shown in red. 
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Figure 26. 




Yearday 2002 


Cross-shore bulk energy flux plotted against yearday in 2002 from the 
direetional wave buoys offshore of Marina and Sand City, Monterey Bay, USA. 
Mean values are shown in red. Onshore values are negative. 
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Figure 27. Alongshore bulk energy flux plotted against yearday in 2002 from the 
directional wave buoys offshore of Marina and Sand City, Monterey Bay, USA. 
Mean values are shown in red. Southward values are positive, indicating fluxes 
toward Del Monte Beach. 
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The Coastal Data Information Program 
@ Scripps Institution of Oceanography 

CA Dept, of Boating and Waterways - U.S. Army Corps of Engineers 
U.S. Naval Postgraduate School 

*** NOTE: WAVE HEIGHT COLOR SCALE CHANGES WITH CONDITIONS *** 
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Figure 28. Output from O’Reilly’s pure refraction wave model for the Monterey Bay. 

The model is initialized with directional spectra shown in the sub-graphic, 
indicating the directional distribution of the energy at various frequencies incident 
at NOAA buoy 46042. Significant wave height (feet) is color coded below the 
graph. 
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Figure 29. 


Figure 30. 


Modeled vs. Actual Significant Wave Height for Sand City. Monterey Bay. USA. 2002 



Yearday 2002 


A eomparison of the significant wave heights at Sand City in 2002, as 
generated by O’Reilly’s refraction model and a directional wave-rider buoy 
located offshore in 18m of depth. 


Modeled vs. Actual Significant Wave Height for Marina. Monterey Bay. USA. 2002 



Yearday 2002 


A comparison of the significant wave heights at Marina in 2002, as generated 
by O’Reilly’s refraction model and a directional wave-rider buoy located offshore 
in 18m of depth. 
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Significant Wave Height (m) Regression for Sand City, Monterey Bay, USA, 2002 



Figure 31. A linear regression of the reeorded by the Sand City wave buoy with 
O’Reilly’s refraction model run at the same location and depth. 


Significant Wave Height (m) Regression for Marina, Monterey Bay, USA, 2002 



Figure 32. A linear regression of the recorded by the Marina wave buoy with 
O’Reilly’s refraction model run at the same location and depth. 
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Figure 33. Alongshore wavelength, L^, plotted against the direetional spreading 
parameter, , as calcuated by the Delft 3-D morphodynamie model (Renters et 
al. 2003). Model runs are for constant = \m and 6^^^^ = 0° and an initially 
longshore uniform barred beach. At = 0°, the circulation is self-organized 
and a minimum L^. is predicted. Adding a small amount of directional spreading, 
(7bulk = 2°, made the beach quasi-forced and resulted in a maximum L^,. Further 
increases in a bulk resulted in a decrease in . 
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